HOMOGENIZED MODELS FOR A SHORT-TIME FILTRATION AND 
FOR ACOUSTIC WAVES PROPAGATION IN A POROUS MEDIA 

ANVARBEK MEIRMANOV* 

Abstract. We consider a linear system of differential equations describing a joint motion of elas- 
tic porous body and fluid occupying porous space. The rigorous justification, under various conditions 
imposed on physical parameters, is fulfilled for homogenization procedures as the dimensionless size 
of the pores tends to zero, while the porous body is geometrically periodic and a characteristic time 
of processes is small enough. Such kind of models may describe, for example, hydraulic fracturing 
or acoustic or seismic waves propagation. As the results, we derive different types of homogenized 
equations involving non-isotropic Stokes system for fluid velocity coupled with acoustic equations for 
the solid component or different types of acoustic equations, depending on ratios between physical 
parameters. The proofs are based on Nguetseng's two-scale convergence method of homogenization 
in periodic structures. 
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1. Introduction. The paper addresses the problem of a joint motion of a de- 
formablc solid {elastic skeleton), perforated by system of channels or pores {porous 
space) and a fluid, occupying porous space. In dimensionless variables (without 
primes) 

x' ^ Lx, t'^Tt, w'^Lw, p'^^pqPs, p'j:=poPf, F'^gF, 

differential equations of the problem in a domain g R'^ for the dimensionless dis- 
placement vector w of the continuum medium have a form: 

(1.1) arP-g-^=divxP + pF, 

(1.2) P = xa^,D[x, — j + (1 - x)(yxD{x, w)-{q + tt)/, 



(1.3) g=p- 



oty dp 
ap dt ' 



(1.4) p + xapdiv^w = 0, 

(1.5) vr + (1 - x)a,,diva;U; = 0. 

The problem is endowed with homogeneous initial and boundary conditions 

(1.6) -u;|i=o = 0, ^|t=o = 0, xen 



(1.7) w^O, xeS = dn, t>0. 

Here and further we use notations 

D{x, u) = (1/2) {V,u + (V,w)^) , p^xpf + i^- x)Ps. 
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In this model the characteristic function of the porous space x(a;) and a dimensionless 
vector F{x,t) of distributed mass forces are known functions. 

Dimensionless constants {i — t, v, . . .) are defined by the formulas 

L h' 2fi 



gr^ TLgpo rLgpo 

_ c^pf _ V _ 2X 

Lg LgpQ Lgpo 

where p is the viscosity of fluid or gas, is the bulk viscosity of fluid or gas, A and rj 
are elastic Lame's constants, c is a speed of sound in fluid, L is a characteristic size 
of the domain in consideration and r is a characteristic time of the process. 
For more details about Eqs. (frT|) - (fTS]) see [7], [2], [TO] . 

Our aim is to derive all possible limiting regimes (homogenized equations) for the 
problem pTT]) - pTT)) as e \ 0. 

To do that we accept the following constraints 

Assumption 1.1. domain fl — (0, 1)'^ is a periodic repetition of an elementary 
cell — eY , where Y = (0, 1)'^ and quantity 1/e is integer, so that fl always contains 
an integer number of elementary cells Yf . Let Yg be a "solid part" of Y , and the 
"liquid part" Yf - is its open complement. We denote as ^ — dYf D dY^ and 7 is 

-surface. A porous space Sl^ is the periodic repetition of the elementary cell sYf, 
and solid skeleton fll is the periodic repetition of the elementary cell sYs . A boundary 

~ dQl n cJSly: is the periodic repetition in Q of the boundary e^. The "solid 
skeleton" £7^ and the "porous space" 17^ are connected domains. In these 
assumptions 



P = P'{x) - x'{x)pf + (1 - X'{x))ps. 

Let e be a characteristic size of pores I divided by the characteristic size L of the 
entire porous body: 

I 

Suppose that all dimensionless parameters depend on the small parameter e and 
there exist limits (finite or infinite) 

lima^(£) = ^o, hmQ;A(e) = Ao, lim Q;^(e) = tq, 

eSO e\0 £\,0 

lim Q;^(e) = 770, lim ap{e) = p*, lim a^{e) — vo, 

e\0 £\,0 £\,0 

lim ^ — pi lim ^ = Ai. 

e\o £\o 

The first research with the aim of finding limiting regimes in the case when the 
skeleton was assumed to be an absolutely rigid body was carried out by E. Sanchez- 
Palencia and L. Tartar. E. Sanchez-Palencia [TUl Sec. 7.2] formally obtained Darcy's 
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law of filtration using the method of two-scale asymptotic expansions, and L. Tar- 
tar [ini Appendix] mathematically rigorously justified the homogenization procedure. 
Using the same method of two-scale expansions J. Keller and R. Burridge [5] derived 
formally the system of Biot's equations from the problem (|l.ip - (|1.7p in the case when 
the parameter was of order e^, and the rest of the coefficients were fixed indepen- 
dent of e. Under the same assumptions as in the article [2], the rigorous justification 
of Biot's model was given by G. Nguetseng [9] and later by Th. Clopeaut et al. [3]. 
The most general case of the problem (|l.f p - (|1.7p when 

has been studied in [7]. 

All these authors have used Nguetseng's two-scale convergence method [HIIS]- 
In the present work by means of the same method we investigate the rest of all 
possible limiting regimes in the problem (jl.ip - (|1.7p . Namely, if tq = oo, which is a 
case of short-time processes, then re-normalizing the displacement vector by setting 

we reduce the problem to the case tq = 1 and /io < oo. Here the only one case Aq = 
needs an additional consideration. 

Therefore we restrict ourself by the case, when 

j/Q, Mo < oo; Ao = 0, tq = I, < p,, 77Q. 

We show that in the case /io > the homogenized equations are non-isotropic 
Stokes equations for fluid velocity coupled with acoustic equations for the solid com- 
ponent, or non-isotropic Stokes system for the one- velocity continuum (theorem l2.2p . 
In the case /xq = the homogenized equations are different types of acoustic equations 
for two- velocity or one- velocity continuum (theorem 12. 3p . 

2. Main results. As usual, equation (|f .ip is understood in the sense of distri- 
butions. It involves the proper equation (II. ip in a usual sense in the domains 17^ and 
ri^ and the boundary conditions 

(2.1) [w] =0, x^e t > 0, 

(2.2) [P] =0, xoeT", t>0 

on the boundary F^, where 

[(p]{xo) = 'f{s){xo) - ip{f){xn), 
f{s){xo) = lim ip{x), ifi(j){xa) = lim ifi{x). 

=0 e Si; n, e n'j, 

There are various equivalent in the sense of distributions forms of representation 
of equation (|I.ip and boundary conditions (|2.ip - p.2p . In what follows, it is convenient 
to write them in the form of the integral equalities. 

We say that functions {w^,p^, , tt^) are called a generalized solution of the prob- 
lem (|l.ip - (|1.7p , if they satisfy the regularity conditions 



(2.3) 
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in the domain Qt — x (0, T), boundary conditions p.7p in the trace sense, equations 
p.3p - (|1.5p a.e. in and integral identity 



{(1 - x')axD{x, w^) - [q" + vr^)/} : ^(a;, if)] dxdt ^ 



(2.4) 



for all smooth vector- functions if = (p{x,t) such that (p\dn = 'p\t=T = d(p/dt\t=T = 0. 

In (|2.4p by A : i? we denote the convolution (or, equivalently, the inner tensor 
product) of two second-rank tensors along the both indexes, i.e., A : B — tr {B* oA) = 

In what follows all parameters may take all permitted values. If, for example, 
tjq^ — , then all terms in final equations containing this parameter disappear. 



The following theorems 1 2 . 1H2 . 31 are the main results of the paper. 

Theorem 2.1. Let F, dF/dt and d'^F/dt^ are hounded in L'^{n). Then for all 
e > on the arbitrary time interval [0, T] there exists a unique generalized solution of 
the problem 11. |i.7| ) and 

(2.5) ,f,fX-^Whn<Co, 



(2.6) max^||^;^-|V, — (t)| + (l-x^)V^|V. — (i)|||2,o < Co, 



(2.7) Wq'h.n^ + Wfh.n^ + — ||^||2,o. < Co 

ap ot 

where Cq does not depend on the small parameter e. 

Theorem 2.2. Assume that the hypotheses in theorem \2.1\ hold, and fiQ > 0. 
Then functions dw^ /dt admit an extension from 51^ = 17^- x (0,T) into 57t such 
that the sequence {t'^} converges strongly in L'^{Q,t) and weakly in L^((0,T); M^2^(r2)) 
to the function v. At the same time, sequences {w^}, {(1 — x)"^^}; {P^}; {(f}; ^'^^ 
{tt'^I converge weakly in L'^(Qt) to w, w'^ , p, q, and tt, respectively. 

I) If Xi = oo, then dw'^ /dt = (1 — m)v = (1 — m)dw/dt and weak and strong 
limits q, p, tt and v satisfy in fix the initial-boundary value problem 

(2 8) ^~di^ diVx{poA{, : D{x,v) + B^TT + BldiVxV+ 

lo - T)diVxv{x, T)dT} ~ V{q + tt) + pF, 



(2.9) 



P* ^dp/dt + Cq : D{x, v) + Cq-tt + (a{ + m)diVxV 

+ /g ^{{t ~ T)diVxV{x,T)dT = 0, 



1^0 dp 1 dp 1 dn , 
2.10 q = P+—^, —^ + -j- + diVxV = Q, 

dt p^ dt rjo dt 

where p = mpf + (1 — m)ps, ni — Jy xdy and the symmetric strictly positively defined 
constant fourth-rank tensor Aq, matrices Cq , Bq , B( and -B2 (^) ^'^'^ scalar s a^, a{ 
and a2{t) are defined below by formulas i5.29]) . iS.Sl]) - iS.SS^) . 
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Differential equations i2.8]) are endowed with homogeneous initial and boundary 
conditions 

(2.11) t;(a;,0) = 0, x e il, v{x,t) = 0, x e S, t > 0. 

II) If Xi < oo, then weak and strong limits , q, p, tt and v satisfy in VLt the 
initial-boundary value problem, which consists of Stokes like system 



(2.12) 



(2.13) 



Pfmdv/dt + psd'^w'^ /dt^ + V((j' + tt) — pF = divxiB^-K^ 
PqAq : D{x, v) + B(diVxV + J* B^it — T)divxv{x, T)dT}, 

P*^dp/dt + Cq : D{x, v) + a^-n: + (a{ + m)diVxV "1 
+ a{(t — T)divxv{x, T)dT — 0, J 



(2.14) '^-P+'t 

P* at 

for the liquid component coupled with a continuity equation 

1 dp 1 dn dw'^ 

(2.15) — -TTT + divx—- + mdm^v ^Q, 

at rjo ot at 

the relation 

dw' /■* 

(2.16) -J— ^ {l-m)v{xA) + B^it-T) ■ z{x,T)dT, 

at Jo 

1 du 

Z{x,t) = -~- VxTT{x,t) + PsF{x,t) - ps — {x,t) 

1 — m ot 
in the case of Ai > 0, or the balance of momentum equation in the form 
d^w^ dv 1 

(2.17) ps-^ - psB^, • + " - • ("T^^-'' + P^^^ 

in the case of Xi — for the solid component. The problem is supplemented by 
boundary and initial conditions \2.11\) for the velocity v of the liquid component and 
by the homogeneous initial conditions and the boundary condition 

(2.18) w''{x,t) -nix) ^Q, {x,t) e S, t > 0, 

for the displacement w'^ of the solid component. In Eqs. \2.16\) - l2.18\) n[x) is the 
unit normal vector to S at a point x € S, and matrices B\{t) and i?| are given 
below by Eqs. i f 5. 39) and \5.4^1^ , where the matrix ((1 — rn)/— i?|) is symmetric and 
positively definite. 

Theorem 2.3. Assume that the hypotheses in Theorem \2.1\ hold, and 



Po = 0; p*, 770 < 00. 
Then there exist functions 11;^, w^^ € L°°(0, T; ^^2^(17)) such that 

w'j = in X (0, T), wl = w' m x (0, T) 
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and sequences {p^}, {q^}, {tt"^}, {w^}, {x'^'^'^}j i'^^f} ^^'^ {''^sl converge weakly in 
L'^(JIt) to functions p, q, tt, w, , w f and Wg respectively as £ \ 0. 

I) If Hi = Xi = oo, then Wf ^ Wg = w and functions w, p, q and tt satisfy in 
Ot the system of acoustic equations 

d^w 1 

(2.19) ^at^--(T3^^^- + ^^' 

(2.20) —p + — TT + diVa:W = 0, 

P* Vo 

(2.21) <i = p + -% -^-r^-> 

p^ Ot m 1 — m 

homogeneous initial conditions 

dw 

(2.22) w{x, 0) = — (a;, 0) = 0, xeil 
and homogeneous boundary condition 

(2.23) w{x, t) ■ n{x) = 0, x e S,t>0. 

II) If Hi = oo and Xi < oo, then functions Wf — w, , p, q and tt satisfy in fix the 
system of acoustic equations, which consist of the state equations 112. 21\) and balance 
of momentum equation 

(2.24) p,^^ ~ (T^^^^ + 
for the liquid component, continuity equation 

(2.25) — pH TT + mdivxWf + divxW^ = 0, 

P* Vo 

and the relation 

(2.26) ^ = (1-^)^ + I^Btit-r)- z%x, T)dr, 

z''{x,t) = ---^—VxT^{x,t)+pgF{x,t) - ps^-^{x,t) 
1 — m ot'' 

in the case of Xi > 0, or the balance of momentum equation in the form 

(2.27) Ps^^ PsB^2 ■ ^ + ((1 - m)I SI) • (- Y^V.TT + p,F) 

in the case of Xi = for the solid component. The problem h2.21\) . {2.24^ - 12.27]) is 
supplemented by homogeneous initial conditions h2.22\) for the displacements in the 
liquid and the solid components and homogeneous boundary condition h2. 2S\] for the 
displacements w = mw f + to* . 
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In Eas A2.2^) - \2.2l^ matrices B1{t) and B?2 are the same as in theorem[ 

III) // /ii < oo and Ai = oo , then functions , — w, p, q and tt satisfy in 
Qt the system of acoustic equations, which consist of the state equations V2.21]) and 
the balance of momentum equation 

d^w^ d'^w 1 

(2.28) + - H V - - (1^^^^ + 

for the solid component, the continuity equation 

(2.29) — p H vr + diVxW^ + (1 — m)diVxWs — 0, 

P* Vo 

and the relation 

(2.30) ^ = + B({t - r) ■ zf{x, r)dr, 

zf{x,t) = -—Vxq{x,t) + p}F{x,t) - pf-^{x,t) 
in the case of Ai > 0, or the balance of momentum equation in the form 

(2-31) P/^ = PfB', ■^ + iml- Bi) ■ (-Iv.g + pfF) 

in the case of Xi =0 for the liquid component. The problem i2.21\) . i2.28\) - ^2.31\) 
is supplemented by homogeneous initial conditions i2.22\) for the displacements in the 
liquid and the solid components and homogeneous boundary condition h2. 23\) for the 
displacements w = + (1 — m)Wg. 

In Eqs.lMB-lMW matrices B{{t) and B{ are given below by formulas \6. 32\) - 
16'. 33]) . where the matrix {ml — B^) is symmetric and positively definite. 

IV) If pi < oo and Ai < oo, then functions w, p, q and tt satisfy in the 
system of acoustic equations, which consist of the continuity and the state equations 
i2.20\) and h2.21]} and the relation 



(2.32) j\-^t~T)- Virix, T)dT + f{x, t), 



where B'^ it) and f{x,t) are given below by Eqs. ( \67^^ and {6.^1^ . 

The problem \2.2U\) . h2.21\) . i2.32^) is supplemented by homogeneous initial and 
boundary conditions V2.22\) and V2.23\) . 

3. Preliminaries. 

3.1. Two-scale convergence. Justification of theorems I2.1H2.31 relies on sys- 
tematic use of the method of two-scale convergence, which had been proposed by 
G. Nguetseng [S] and has been applied recently to a wide range of homogenization 
problems (see, for example, the survey [S|). 

Definition 3.1. A sequence {<p'^} C L^{Q,t) is said to be two-scale convergent to 
a limit Lp £ IF'{Q,t x Y) if and only if for any 1-periodic in y function a — cf(x, t, y) 
the limiting relation 



(3.1) 



lim / tp'^{x,t)a {x,t,x/e) dxdt ^ / / ip{x,t,y)a{x,t,y)dydxdt 
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holds. 

Existence and main properties of weakly convergent sequences are established by 
the following fundamental theorem 'W, '6' : 

Theorem 3.2. (Nguetseng's theorem) 

1. Any bounded in L^{Q) sequence contains a subsequence, two-scale convergent 
to some limit ip G L^(r2y x Y). 

2. Let sequences {(^"^j and {eVj-^s^} be uniformly bounded in L'^{Q.t)- Then there exist 
a 1-periodic in y function ip — ip{x,t,y) and a subsequence {9?^} such that ip,'S/yip G 
L^{Q,T X Y), and pf and eS/ xf^ two-scale converge to (p and "S/yip, respectively. 

3. Let sequences {(p^} and {Vj^^s*^} be bounded in L^{Q). Then there exist functions 
ip e L^(f2'r) and ip € L^ijlx x Y) and a subsequence from {(p^^ such that ijj is 1- 
periodic in y, Vyijj G L^{Q,t x Y), and 95^ and "S/ x'f^ two-scale converge to ip and 
\I xV'ix, t) + \I yi]j{x, t, y), respectively. 

Corollary 3.3. Let a S L'^iY) and cr^{x) :— a{x/e). Assume that a sequence 
{v?^} C L^{rtT) two-scale converges to ip E L^{^It x Y). Then the sequence u^ip^ 
two- scale converges to aip. 

3.2. An extension lemma. The typical difficulty in homogenization problems 
while passing to a limit in Model i?^ as e \ arises because of the fact that the 
bounds on the gradient of displacement 'SI may be distinct in liquid and rigid 
phases. The classical approach in overcoming this difficulty consists of constructing 
of extension to the whole of the displacement field defined merely on fig. The 
following lemma is valid due to the well-known results from [1] |4] . We formulate it in 
appropriate for us form: 

Lemma 3.4. Suppose that assumvtion ll .l\ on geometry of periodic structure holds, 
€ W2i^l) and = on Sg = dil^ fl dft in the trace sense. Then there exists 
a function G W2^(ri) such that its restriction on the sub-domain fl^ coincide with 
■0% i.e., 

(3.2) {l-x'ix)){a'ix)^i;%x))^0, x e n, 

and, moreover, the estimate 

(3-3) ||a^||2,o < C\\r\\2,n^^, HV^a^l^.n < C\\\/xr\kni 

hold true, where the constant C depends only on geometry Y and does not depend on 
e. 

3.3. Friedrichs Poincare's inequality in periodic structure. The follow- 
ing lemma was proved by L. Tartar in 10, Appendix] . It specifics Friedrichs-Poincare's 
inequality for e-periodic structure. 

Lemma 3.5. Suppose that assumptions on the geometry of hold true. Then 

o 

for any function ip & W2 {^f) the inequality 

(3.4) [ \ipfdx < Ce^ [ IVx'pfdx 

holds true with some constant C , independent of e. 
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3.4. Some notation. Further we denote 
1) 



JY JY JY 



{(p)n= / (pdx, {ip)nT = / ipdxdt. 
Jn jut 

2) If a and b are two vectors then the matrix a (X) b is defined by the formula 

(a (g) b) • c = a(b • c) 

for any vector c. 

3) If B and C are two matrices, then B (g) C is a forth-rank tensor such that its 
convolution with any matrix A is defined by the formula 

{B(g)C):A = B{C : A) 

4) By we denote the 3 x 3-matrix with just one non-vanishing entry, which is 
equal to one and stands in the i-th row and the j-th column. 

5) We also introduce 

^ _^ = l(e, ® + Bj ® ei), 

where (ei, 62, 63) are the standard Cartesian basis vectors. 

4. Proof of theorem [2TT1 Estimates (P3 |) - ((^ follow from 

^mM^(^/a;;||diV;:,5i(;''/ai(t)||2,n5 + ^/ooiW^ xdw^ / dt{t)\\2^Q.'=^ 



(4.1) +^Ww^lde{t)\\2^n + ya^||div,9«;Vat(t)||2,n=) 



where Cq is independent of e. Last estimates we obtain if we differentiate equation for 
with respect to time, multiply by d^w^ /dt^ and integrate by parts using continuity 

and state equations (|1.3p - (|1.5p . The same estimates guarantee the existence and 

uniqueness of the generalized solution for the problem (jl.ip - p.7p . 

li p.t + rjo < 00, then estimate ()2.7p for pressures follows from estimate (|4.H) and 

continuity and state equations (|1.3p - (jl.Sp . 

For the case -I- 770 = 00 estimate (|2.7p follows from integral identity (|2.4p and 

estimates (I4.ip as an estimate of the corresponding functional, if we re-normalized 

pressures, such that 



/ {q^x,t) + T:%x,t))dx = 0. 
Jn 

Indeed, integral identity (|2.4p and estimates (|4.ip imply 

I / iq' +^r')diyx'>pdx\<C\\V^P\\2,n■ 
Jn 
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Choosing now ■0 such that {q^ + tt^) = div^ip we get the desired estimate for the 
sum of pressures {q^ + tt^). Such a choice is always possible (see [S]), if we put 

ip ^\/(p + ipo, divjjV'o = 0, A(/7 = (j^ + tt", (p\dn = 0, C^f + il>o)\an^O- 

Note that the re-normalization of the pressures (q^ +7r^) transforms continuity equa- 
tions (|1.4p - (|1.5|) for pressures into 

(4.2) —p' + x'diY^w' = -fi'x', 

Up m 

(4.3) — TT- + (1 - x")div,t.;^ = -^—pe^i _ ^e^^ 

(1 — m) 

where 

= (x"div,u;")n. 

In what follows we will use equations (|4.2p and (|4.3p only if -I- 770 = cxd . 

Note that for the last case the basic integral identity (|2.4p permits to bound only 
the sum (g*^ -I-tt^). But thanks to the property that the product of these two functions 
is equal to zero, it is enough to get bounds for each of these functions. The pressure 
is bounded from the state equation (|1.3p . if we substitute the term {a,^ / ap)dp^ / dt 
from the continuity equation 14.21 and use estimate (14. ip . 

5. Proof of theorem Hm 

5.1. Weak and two-scale limits of sequences of displacement and pres- 
sures. On the strength of theorem 12. 1[ the sequences {p^}, {q'^}, {tt^} and {w'^} 
are uniformly in e bounded in L'^^flx)- Hence there exist a subsequence of small 
parameters {s > 0} and functions p, q, tt and w such that 

(5.1) p^ P, <f ^ 1, TT^ TT, ^ W 

weakly in L^(Ot) as e \ 0. 

Moreover, due to lemma [3T4l there is a function v'^ e L°°{0, T; W2(^)) such that 
t)^ — dw^/dt in ri/ X (0,T), and the family {v'^} is uniformly in e bounded in 
L°°(0, T; T4^2^(ri)). Therefore it is possible to extract a subsequence of {e > 0} such 
that 

(5.2) ^vwciMyinL^{Q,T;Wl{VL)) 

as e \ 0. 

Note also, that 

(5.3) {\-x')axD{x,w')^Q. 

strongly in L'^{Q,t) as e \ 0. 

Relabeling if necessary, we assume that the sequences converge themselves. 

On the strength of Nguetseng's theorem, there exist 1-periodic in y functions 
P{x,t,y), Il{x,t,y), Q{x,t,y), W{x,t,y) and V{x,t,y) such that the sequences 
{p^}, {tt^}, {g"}, {if^} and {Va;t>^} two-scale converge to P{x,t, y), n(a;, t,y), Q{x, t,y), 
W{x, t, y) and ^ xV + Vy V(a;, i, y), respectively. 

o 

Note that the sequence {dYv^w^} weakly converges to div^,!!; and v E L^{0, T; W2 
(fi)). Last assertion follows from the Friedrichs-Poincare's inequality for v'^ in the e- 
layer of the boundary S and from convergence of sequence {v^} to v strongly in 
L'^inr) and weakly in L^{{0,T);W^{n)). 
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5.2. Micro- and macroscopic equations I. We start this section with macro- 
and microscopic equations, connected with continuity equations. 

Lemma 5.1. For all x € Q and y (z Y weak and two-scale limits of the sequences 
{p^}, {tt^}, {q^}, {w^}, and {v"^} satisfy the relations 



(5.4) n-7r(l-x)/(l-m); 

(5.5) q=p + iyop-^dp/dt, Q = P + lyop-'dP/dt; 

(5.6) p^^dp/dt + mdiVxV + {diVyV)Yf — dp/dt; 

(5.7) p-^dP/dt + xidivxV + diVyV) = {x/m)d(3/dt; 

(5.8) p/p* + n/rjo + divxW = 0; 

(5.9) w{x, t) ■ n{x) ^0, xe S,t>0; 

(5.10) diVyW = 0; 

(5.11) dW/dt = xv + {l-x)dW/dt, 



where (3 = {{diVyU)Yf)n, if P* + ?7o = cxd and P = 0, if p^ + 770 < oo and n{x) is the 
unit normal vector to S at a point x € S . 

Proof. In order to prove Eq. (|5.4p . into Eq. (|2.4p insert a test function xp"^ ~ 
exp {x,t,x/e), where xp(x,t,y) is an arbitrary 1-periodic and finite on Yg function in 
y. Passing to the hmit as e \ 0, we get 

(5.12) V;;n(a;,t,y) = 0, y G F,. 

Next, fulfilhng the two-scale limiting passage in equality 

x'tt' = 

we arrive at 



xn = o 

which along with Eqs. (|5.12p justifies Eq. (|5.4p . 

Eqs. (|5.5p - (|5.9p appear as the results of two-scale limiting passages in Eqs. (|1.3p . 
(|4.2p - (|4.3p with the proper test functions being involved. Thus, for example, Eqs. 
(|5.8p and (|5.9p arise, if we consider the sum of Ea. (|4.2p and Eq. (|4.3p 

(5.13) —p' + —tt' + diY^w' = — i -13'ix' - m), 

Up m[l — to) 

multiply by an arbitrary function, independent of the "fast" variable x/e, and then 
pass to the limit as e \ 0. In order to prove Eq. (|5.10p . it is sufficient to con- 
sider the two-scale limiting relations in Eq. (|5.13p as e \ with the test functions 
£"0 i^/s) h(x, t), where %p and h are arbitrary smooth test functions. In order to prove 
Eq. (|5.1ip it is sufficient to consider the two-scale limiting relations in 

X^'idw^ /dt-v^) = 0. 

□ 

Corollary 5.2. If p^ + rjo — 00 , then weak limits p, tt and q satisfy relations 



(5.14) 



{p)n = (7r)o = {q)n = 0. 
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Lemma 5.3. For all {x,t) G fiy the relations 

(5.15) d^Vy{^lox{Diy, V) + D{x, v)) - (Q + ^^^^ • /} - 0, 

(1 — m) 

holds true. 

Proof. Substituting a test function of the form ■0^ = exp {x,t,x/e), where 
xp{x,t,y) is an arbitrary 1-periodic in y function vanishing on the boundary 5, into 
integral identity (|2.4p . and passing to the hmit as e \ 0, we arrive at Eq. (|5.15p . □ 

Lemma 5.4. Let p — mpj + (1 — m)ps. Then functions ~ (W^)"^; 1 "■''T-d 
TT satisfy in fix the system of macroscopic equations 

dv d^w^ 

(5.16) pfm— + ps-^ - pF = 

div^{po{mD{x, v) + {D{y, V))y^) - {q + n) ■ I}. 



Proof. Eqs. (|5.16p arise as the hmit of Eqs. ()2.4p with test functions being finite 
in Qt and independent of e. □ 

5.3. Micro- and macroscopic equations II. In this section we derive macro- 
and microscopic equations for the sohd component. 

Lemma 5.5. If Xi = oo, then the weak limits of {v^} and {dw'^/dt} coincide. 

Proof. Let "9{x,t,y) be an arbitrary function periodic in y. The sequence {f3^}, 
where 



/ \/axS/w^ {x,t)'^{x,t,x/ £)dx, 

is uniformly bounded in e. Therefore, 

/ eVu;''*(a;,t,a;/e)dx = ^=/3= 

as e \ 0, which is equivalent to 

/ / W{x,t,y)\Iy^{x,t,y)dxdy ^0, 
Jn Jy 

or W{x, t, y) = w{x, t). □ 

Lemma 5.6. Let Ai < oo. Then the weak and two-scale limits tt and W satisfy 
the microscopic relations 



d'^W 1 

(5.17) Ps^-^^Xil^yW -VyR-- V.,TT + psF, yeYs, 

ot^ 1 — m 



dW 

(5.18) "df^'"' ^ ^ ^ 

in the case Ai > 0, and relations 



d^W 1 
(5.19) p,— = -Vyi?- S/^TT + psF, yeY, 

Ot-^ 1 — TO 
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(5.20) (__-„).n = 0, ye J 

in the case Ai = 0. 

In Eq. \5.20\) n is the unit normal to 7. 

Proof. Differential equations (|5.17p and (|5.19p follow as e \ from integral 
equality (|2.4[) with the test function xjj = (p{xe^^) ■ h{x, t), where ip is solenoidal and 
finite in K,.. 

Boundary conditions (|5.18p are the consequences of the two-scale convergence of 
1 1 
{a^V xW^} to the function WyW{x, t, y). On the strength of this convergence, the 

function V yW{x,t,y) is i^-integrablc in Y. The boundary conditions (|5.20p follow 

from Eqs. (|5TU)) - ([5ll|) . □ 

5.4. Homogenized equations I. Here we derive homogenized equations for 
the liquid component. 

Lemma 5.7. // Ai =00 then dw/dt = v and the weak limits v, p, q, and tt 
satisfy in VLt the initial-boundary value problem 

dv f f f 

(5 21) ^Ih " divxino^o • D{x,v) + B(^tt + B{diVxV+ 

lo ^2(t - T)div^v{x,T)dT} - V{q + tt) + pF, 



(5.22) 



P* ^dp/dt + Cq : D{x, v) + AqYt + (a{ + m)diVxV 
+ /p 02 (i — T)diVxv{x, T)dT — 0, 



/K r,Q^ , ^odp I dp I d-K 

p* at at rio ot 

where the symmetric strictly positively defined constant fourth-rank tensor A^, matri- 
ces C^,B^, B( and i?2(^) '^'^'^ scalar s a^, a( and 02 (t) are defined below by formulas 

lEM, ISM - WM- 

Differential equations 115. 21\) are endowed with homogeneous initial and boundary 
conditions 

(5.24) v{x,0) = 0, xen, v(x,t) = 0, xeS, t>0. 

Proof. In the first place let us notice that v = dw/dt due to lemma 
The homogenized equations (|5.2ip follow from the macroscopic equations (|5.16p . 
after we insert in them the expression 

PQ{D(y, V))Yf = HoM : D(x, v) + B^n + B(divxV + J B({t- T)div^v{x, T)dT. 

In turn, this expression follows by virtue of solutions of Eq. (j5.5p in the form 

Q = P - i'ox(divxi' + divyV) + iy(){x/m)df3/dt 
and Eas. (|5.6p and (|5.15p on the pattern cell Yf. Indeed, setting 
3 
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+ / V^^\y,t-T)div^v{x,T)dT, 
Jo 



3 



+ / <3^^^ (y. * - T)divxw(x, r)dr, 
Jo 

Jo 



where 



we arrive at the following periodic-boundary value problems in Y: 

divyixiDiy, + - Q(^^) . /)} = 0, 1 

(5.25) noQ^'J^ + PoXdivyV^'^'^ = for < oo, 
Xdivy y^^-''^ = for = oo; 

div,{Mox^(y, - (0(°) + • /} = 0, 

(5.26) Q(0) + j.o;^divj,F(°) = for p» < oo, 
Xdivy V^°^ =0 for = oo; 



div,{Mox£'(y,vW)-Q^'^-/} = o, 

(5.27) + i/ox(divyF(^) + 1) = for < oo, 
x(diVyF(^^ + 1) = for p^ = oo; 

^ 0, Q(2) = p(2) = for = oo and 

divjMox^(y,'^^^'^)-Q^'^-/} = o, 

(5.28) Q(2) ^p(2) _^^^divj,y(2), 
lMap(2)/5i + xdivj,V(2) =0, P(2)(y,0) =p. ^ 

for < oo. 

Note that for + % = oo 

3 

/3 = ((div,V)a)y, = ^ (div,V(*^))y,(A,)n+ 
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(div^y(°V, Wo + (div^y(iV/(div.t^>o = 

due to homogeneous boundary conditions for v{x,t) and Eq. (|5.14p . 

On the strength of the assumptions on the geometry of the pattern "Hquid" ceU 
Yf, problems (|5.25p - (|5.27p have unique solution, up to an arbitrary constant vector. 
In order to discard the arbitrary constant vectors we demand 

Thus 

3 3 

(5.29) 4 = J"' J"' + 4^Y1 ® J"'- 

Symmetry of the tensor Aq follows from symmetry of the tensor A(. And symmetry 
of the latter one follows from the equality 

(5.30) {D{y, y(*^")))y, : = -{D{y, V^^'^) : D{y, y('=')))y, 

-sgn( — )— divyy('J)divj,y(^'\ 
P* Mo 

which appears by means of multiplication of Eq. (lOSj) for V'-'^'^ by F^'''^ and by 
integration by parts using the continuity equation. 

This equality also implies positive definiteness of the tensor Aq. Indeed, let C be 
an arbitrary symmetric matrix. Setting 

and taking into account Eq. (l5.30p we get 

{Diy,Z))Y, : C = -{D{y,Z) : i?(y, Z))y, - sgn(l)^(div,Z)2, 

P* Po 

This equality and the definition of the tensor Aq give 

(4 : C) : C = my, Z) + C) : {D{y, Z) + C))n + sgn(l)^(div,Z)^ 

P* Aio 

Now the strict positive definiteness of the tensor A^ follows from the equality 
immediately above and the geometry of the elementary cell Yf . Namely, suppose that 
(Ag : C) : C = for some function C, such that C : C = 1- Then {D{y, Z) + () = 0, 
which is possible iff Z is a linear function in y. On the other hand, all linear periodic 
functions on Yf are constant. Finally, the normalization condition (V^*-''')y^ = 
yields that Z = 0. However, this is impossible because the functions v'^^-'^ are linearly 
independent. 

Finally, Eqs. ([5?22l) and ([5?23l) for the pressures follow from Eqs. ^M), 
(|5?8l) and 

{divyV)Yf = Cq : D{x,v) + qItt + a(divxV + / a^it — T)divxv{x,T)dT 
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with 

(5.31) b{ ^fio{Diy,V^'^))Yf, * = 0,1,2, 



(5.32) = J2 {diyyV^''%r', a{ = (div^V z = 0, 1, 2. 
□ 

5.5. Homogenized equations II. We complete the proof of theorem 12.21 with 
homogenized equations for the sohd component. 

Let Ai < oo. In the same manner as above, we verify that the hmit v of the 
sequence {f^} satisfies the initial-boundary value problem likes (|5.2ip - (|5.24p . The 
main difference here that, in general, the weak limit dw/dt of the sequence {dw'^ /dt} 
differs from v. More precisely, the following statement is true. 

Lemma 5.8. // Ai < oo then the weak limits v, to", p, q, and tt of the sequences 
{d^}, {(1 — x'^)w'^}, {p^}, {q"^}, and {tt^} satisfy the initial-boundary value problem 
in fix, consisting of the balance of momentum equation 

dv d^w'^ 

(5.33) pfm— + Ps-g-2- + V(9 + n)-pF = 

div^ifio^o '■ ^{x,v) + B^TT + b( div^v} + / B2{t — T)diVxV{x,T)dT}, 







the continuity equation h5. 22\) and first state equation in i5.23\) for the liquid compo- 
nent, where Aq, -Bg - -B| are the same as in h5.21]) . the continuity equation 

/r o.x I dp I dn dw'' 

(5.34) — Trr ^ kt + div^—— + mdiv^v = 0, 

p* ot 7/0 ot at 

the relation 

(5.35) -—- = {l-m)v{x,t)+ I Bl{t-T) ■ z{x,T)dT, 

<jt Jo 

1 di) 

Z{x,t) = ^.j;TT{x,t) + PsF{x,t) - Ps — {x,t) 

1 — TO ot 

in the case of Ai > 0, or the balance of momentum equation in the form 



d^w^ dv 1 

ot'' ot 1 — m 



in the case of Xi =0 for the solid component. The problem is supplemented by 
boundary and initial conditions {5.24-^ for the velocity v of the liquid component and 
by homogeneous initial conditions and the boundary condition 

(5.37) w^ix, t) ■ n{x) ^ 0, {x, t) e S, t > 0, 

for the displacement w" of the solid component. In Eqs. 15. j'5|] - |5.i57[ ) n(a;) is the 
unit normal vector to S at a point x £ S , and matrices Bf (t) and i?| are given below 
by Eqs. ifOPl) and (^Jl^ . 
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Proof. The boundary condition (|5.37p follows from Ea. (|5.9p . equality 



dt dt 



and homogeneous boundary condition for v. 

The same equality and Eq. (|5.8p imply (|5.34p . The homogenized equations of 
balance of momentum (|5.33p derives exactly as before. Therefore we omit the relevant 
proofs now and focus ourself only on derivation of homogenized equation of the balance 
of momentum for the solid displacements w'^ . 

a) If Ai > 0, then the solution of the system of microscopic equations (I5.10|) . 
(I5.17p . and (I5.18|) . provided with the homogeneous initial data, is given by formula 

W= [ {v[x,T)+Bl[y,t-T)- z)[x,T)dT, R = [ Rf{y,t-T)-z{x,T)dT, 
Jo Jo 

in which 



i=l 1=1 

and the functions W^{y,t) and R^{y,t) are defined by virtue of the periodic 
initial-boundary value problem 

(5.38) = 0, y e 7, t > 0; 
W\y,0)^Q, psdWydt{y,0)^e,, y eY,. 

In Eq. (I5.38|) is the standard Cartesian basis vector. 
Therefore 

(5.39) Bl{t) = {^)ySt). 

Note, that equations (|5.38p are understood in the sense of distributions and the 
function Bl{t) has no time derivative at t = 0. 

b) If Ai = then in the process of solving the system (IS.lOp . ()5.19p . and ()5.20p 
we firstly find the pressure R{x,t,y) by virtue of solving the Neumann problem for 
Laplace's equation in Ys in the form 



R{x, t,y) = ^ R,{y)e^ ■ z{x, t), 

i=l 

where R^{y) is the solution of the problem 

(5.40) Ayi?, =0, yeFs; \7yRi ■ n = n ■ Si, y € 7. 

Formula (|5.36p appears as the result of homogenization of Eq. (|5.19p and 

3 

(5.41) S| = ^(Vi?,(y))y, <8)e„ 
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where the matrix ((1 — m)I — _B|) is symmetric and positively definite. In fact, let 
R = X]i=i ^i^i f'^r ^''^y vector ^. Then 

(B • • c = ((c - yRf)Y, > 

due to the same reasons as in lemma 1^771 

Note, that on the strength of the assumptions on the geometry of the pattern 
"solid" cell Ys, problem (|5.38p has unique solution and problem (|5.40p has unique 
solution up to an arbitrary constant. □ 

6. Proof of theorem 12.31 

6.1. Weak and two-scale limits of sequences of displacement and pres- 
sures. Let fiQ = 0. We again use Lemma [3.41 and conclude that there are functions 
w}, 10^ e L°°(0, T; W^{n)) such that 

w^f = in flf X (0, T), wl = in f],. x (0, T). 

On the strength of theorem 12.11 the sequences {p^}, {9^}, {'^'^}, {""^/li 
{y/o^ViOy-}, {w^} and {.y/oAViOj} are uniformly in e bounded in L^(rir)- Hence 
there exist a subsequence of small parameters {e > 0} and functions p, q, tt, w, Wf 
and Wg such that 

(6.1) p"^ — > P, ~* q, ^ TT, ^ w, w^f ^ Wf, wl ^ Ws 

weakly in L?{^It) as e \ 0. 
Note also, that 

(6.2) {l-x'yxD{xMs)^^. x'af.D{x,w})^0 

strongly in L^(17t) as e \ 0. 

Relabeling if necessary, we assume that the sequences converge themselves. 

On the strength of Nguetseng's theorem, there exist 1-periodic in y functions 
P{x,t,y), Il{x,t,y), Q{x,t,y), W{x,t,y), W f{x,t,y) and Ws{x,t,y) such that the 
sequences {p'^}, {tt^}, {q'^}, {W^}, {'^j} and {wl} two-scale converge to P{x,t,y), 
U{x,t,y), Q{x,t,y), W{x,t,y), Wf{x,t,y) and Ws{x,t,y), respectively. 

Finally note, that if /ii = 00 (Ai = 00), then due to Lemma [5.5l the sequence {i"/} 
({ii;^}) converges strongly to wj (wg) and = (VK)y^ mwf { w'^ ~ {W)ys = 
(1 - m)ws). 

6.2. Micro- and macroscopic equations. As before we start the proof of 
theorem with macro- and microscopic equations, connected with continuity equations. 

Lemma 6.1. For all x € n and y &Y weak and two-scale limits of the sequences 
{p^}, {tt^}, {q^}, {if^}, {w^f}, and satisfy the relations 



(6.3) Q = qx/m, P = px/m, H = 7r(l - x)/(l - m); 

(6.4) g/m = 7r/(l - m), q ^ p + vap'^dp/dt; 

(6.5) p/p, + tt/i]o + diVxW = 0; 

(6.6) w{x, t) ■ n{x) = 0, a; e 5, t > 0; 

(6.7) divyW = 0; 

(6.8) W ^xWf + {l-x)W,. 
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Proof. The derivation of Eas. (|6.3p - (|6.8p is the same as the derivation of Eqs. 
(|5.5p . (|5.8p and (|5.10p - (|5.1ip in lemma [5TT1 Thus, for example, the first relation 
in Eq. (|6.4p follow from Eq. (|6.3p and from the strong convergence of the sequence 
{(9^ + tt'^) to + tt), which implies the equality Q + W = q + tt. □ 

Lemma 6.2. For all {x,t) G fir the relation 

(6.9) Pf^+Ps^ = - (T^^^" + 
holds true. 

Proof. Substituting a test function of the form = xp{x,t) into integral identity 
(|2.4p . and passing to the limit as e \ 0, we arrive at Eq. (|6.9p . □ 

Lemma 6.3. Let /ii = 00 and Ai < 00. Then functions {W , Wf, tt} satisfy in Yg 
the system of microscopic equations 

d'^W 1 

(6.10) p^^^ = x,AyW-VyR'-- V,^ + p,F, yeYs, 

ot^ 1 — m 

(6.11) W^Wf, yej 
in the case Ai > 0, and relations 



(6.12) p^^^ = -\/R-- \/.,7r + psF, y€Ys, 

at'^ 1 — m 



(6.13) {W~Wf)-n^O, ye7 

in the case Ai = 0. 

In Eq. i6.13\) n is the unit normal to 7. 

The proof of this lemma repeats the proof of lemma 15.61 

In the same way one can prove 

Lemma 6.4. Let /ii < 00 and Ai = 00. Then functions {W,Ws, q} satisfy in Yf 
the system of microscopic equations 

(6.14) pf——^^,Ay—~WyRf^-W.q + pfF, yeYf, 

(6.15) W = ws, yej 

in the case pi > 0, and relations 

d^W , 1 

(6.16) p^__^_v^i?/__v,q + pyF, y^Yf, 

(6.17) {W -Ws)-n^Q, yej 
in the case /ii = 0. 

Lemma 6.5. Let pi < 00 and Ai < cxd and p = pfX + Ps(l ^ x)- Then functions 
{W, tt} satisfy in Y the .system of microscopic equations 

pd^W/dt^ + 1/(1- m)W^Tr- pF = 1 
^^■^^^ diVy{p,xDiy,dW/dt) + X,{l-x)Diy,W)-RI}. J 
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In the proof of the last lemma we additionally use Nguetseng's theorem, which states 
that the sequence {eD{x, dw^ /dt)} {{eD{x, w^)}) two-scale converges to D{y, dW /dt) 
{D{v,W)). 

6.3. Homogenized equations. Lemmas 16.11 and 16.21 imply 
Lemma 6.6. Let — \i — oo, then Wf = Ws = w and functions w, p, q and vr 
satisfy in fix the system of acoustic equations 

d^w 1 

(6.19) p =- -V,^ + pF, 

OT^ (1 — m) 

(6.20) —p+—TT + diVxW = 0, 

P* Va 

(6.21) q = p-\ — —<}=-, 

P^: Ot m 1 — TO 

homogeneous initial conditions 

dw 

(6.22) w{x, 0) ^ — (a;, 0) ^ 0, x e 
and homogeneous boundary condition 

(6.23) w{x, t) ■ n{x) = 0, x(^S,t>Q. 



Lemma 6.7. Let /ii = oo and Ai < oo. Then functions Wf, id", p, q and tt 
satisfy in Qt the system of acoustic equations, which consist of the state equations 
116.21]} . balance of momentum equation for the liquid component 

(6.24) p,TO^ + = -^Y3^V.. + pF, 
continuity equation 

(6.25) — p -\ TT + mdivxWf + divxw" = 0. 

P* r]o 



and the relation 



(6.26) 



dt 



(1 - to)^ + Bl{t - t) ■ Z%X, T)dT, 



z''{x,t) = ---^—VxT^{x,t)+p,F{x,t) ~ p,^-^(a;,t) 
1 — to ot'' 

in the case of Ai > 0, or the balance of momentum equation for the solid component 
in the form 

(6.27) Ps^^ PsB^ . ^ + ((1 _ m)I B^) ■ (- j^V^vr + p^F) 
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in the case of Ai — 0. The problem i6.21\) . \6.24^ - {6.2'l^ is supplemented by ho- 
mogeneous initial conditions id. 22\) for the displacements in the liquid and the solid 
components and homogeneous boundary condition i6. 23\) for the displacements w = 
mw f + w'^ . 

In Eas. ll6.26\) -{ 6.27\ l matrices Bl{t) and i?| are the same as in theorem \2.2l 



Proof Eq.jO 



} follows directly from Eq. (|6.9p . The continuity equation (|6.25p 



follows from Eq. (|6.5|l if we take into account the equality 



w — mw f + w . 



The derivation of Eqs. (|6.26P " (l6.27p is exactly the same as in lemma [5?8l □ 

Lemma 6.8. Let fj,i < oo and Ai — oo. Then functions , Wg, p, q and tt 
satisfy in ^It the system of acoustic equations, which consist of the state equations 
16.21]) . the balance of momentum equation for the solid component 

D'^njf (p'nj -, 1 

(6.28) + p.(l - m)^ = - (T^^^- + 

the continuity equation 

1 1 f 

(6.29) — p H TT + diVxW-' + (1 — m)diVxWs — 0, 

P* m 

and the relation 

(6.30) ^=m^+ f Bi{t^r)-zf{x,r)dT, 



zf{x,t) ^ ~^Vxq{x,t) + PfF{x,t) - pf^-^{x,t) 

in the case of Ai > 0, or the balance of momentum equation for the liquid component 
in the form 

(6.31) = pfBi .^ + {ml- Hi) ■ (-^V.g + pfF) 

in the case of Ai — 0. The problem h6.21\) . h6.28\) - \6.31\) is supplemented by ho- 



mogeneous initial conditions \6. 22]) for the displacements in the liquid and the solid 
components and homogeneous boundary condition i6. 23\) for the displacements w = 
+ (1 — m)Ws. 

In Eas. l[K^) - l[KJl]) matrices B{(t) and b{ are given below by formulas \6. S2\) - 
Proof. The proof of this lemma repeats proofs of previous lemmas and 

3 

(6.32) B({t)^C£v\y,t))Yf(^e,, 

i=l 



(6.33) i?2^=^(Vi?{(y))y,®e, 
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where functions solve the periodic initial-boundary value problem 

(6.34) P/^-/iiAy' + Vff = 0, div„y* = 0, yeYf,t>0, 

yej,t>0; pfV'{y,0)^e,, y&Yf, 

and functions R{ {y, t) solve the periodic boundary value problem 

(6.35) yeYf, VyRl-n^n-e,, yej. 

Note that as before, the matrix {{ml — B^) is symmetric and positively defined. □ 
The proof of theorem 12.31 is completed by 

Lemma 6.9. Let Let fii < oo and Ai < oo. Then functions w, p, q and tt satisfy 
in Qt the system of acoustic equations, which consist of the continuity and the state 
equations W.2(J\) and \6.21\) and the relation 

(6.36) ^ = j^B^{t-T)- Vtt{x, T)dT + f{x, t), 

where B'" {t) and f{x,t) are given below by Eqs. |6'.^"^ and ^6.41^ - 

The problem 116. 20\) . 116. 21\) . 116. 36\) is supplemented by homogeneous initial and 

boundary conditions h6. 22]) and h6. 2S\) . 
Proof. Let 

^ = /E^^^^^'^- + Wf{y,t-T)F,{x,T)}dT, 



^= I E{^r(y, t - ^)^(^'^) + ^ - T)F,{x,T)}dT, 

where F = T,LlP^^^ and functions {W^V , t) , RUv , t)} and {W( {y , t) , R[ {y , t)} 
are periodic in y solutions of the system 

dWy{fiixD{y,dWl/dt) + X,{l~x)D{y,Wi)~RiL}^ | 
pd^Wl/dt^, diVyWj^O, yeV, t>0, j=n,F, } 

which satisfy the following initial conditions 

dW 1 
(6.38) WUy,0)^0, p^-^{y,Q) = -- e,, xeY, 



dWf 



(6.39) Wf{y, 0)^0, ^^(y,0) = e„ xeY. 

Then the functions W and R solve the system of microscopic equations (|6.7p and 
and by definit ion w = {W)y- Therefore 

(6.40) i3-(i)^^(^(y,<))^^e,, 
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(6.41) {{x,t)^ / Y,{^)y{t--r)F,{x,r)dr. 

Jo at 

The solvability and the uniqueness of problems (|6.37p . (|6.38p or (|6.37p . (|6.39p follow 
directly from the energy identity 

1 BW^ 1 



/■* dW^ dW^- 1 

{f,,D{y, -^) : D{y, -^))y^{r)dr = -;9^ 

for 1 = 1,2,3 and j = tt, _F. 
Here 

p 

As before, equations (|6.37p are understood in the sense of distributions and the func- 
tion B'^ (t) has no time derivative at t = 0. That is why we cannot represent relation 
(|6.36p in the form of the balance of momentum equation, like (|6.19p or (|6.27p . □ 
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